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Abstract. We provide direct proofs of product and coproduct formulae for Schur functions 
where the coefficients (Littlewood-Richardson coefficients) are defined as counting puzzles. 
"q^" . The product formula includes a second alphabet for the Schur functions, allowing in par- 

' ticular to recover formulae of [Molev-Sagan '99] and [Knutson-Tao '03] for factorial Schur 

. functions. The method is based on the quantum integrability of the underlying tiling model. 
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1. Introduction 

Littlewood-Richardson coefficients are important integers related to Schur functions or 
equivalently to the representation theory of the general linear group; they also appear in the 
cohomology of Grafimannians. Interesting combinatorial formulae for them have been given 
in [11, 10]: these coefficients count puzzles, i.e. certain tilings of a triangle where the tiles 
are decorated elementary triangles and rhombi. However all the proofs of this formula are 
fairly indirect e.g. rely on induction. 

From the mathematical physicist's point of view, Littlewood-Richardson coefficients pro- 
vide a challenge. It is well-known [7] that Schur functions are related to (two-dimensional) 
free fermions. However it is not clear how to define Littlewood-Richardson coefficients in 
this framework. In fact the contruction of [10], as will be discussed here, suggests that the 
right way to describes them involves interacting fermions. Interestingly, the physical model 
in question has in fact been studied in the physics literature. As pointed out in [18], it is 
equivalent to a model of random tilings, the square-triangle triangle model, which has been 
the subject of a lot of activity [19, 8, 1]. Note that this equivalence is not particularly useful 
- in order to solve the square-triangle tiling model, one usually goes back to the model of 
tilings of decorated triangles and rhombi. The most important feature of this model for 
our purposes is that it is integrable: the scattering of the elementary degrees of freedom 
(the aforementioned interacting fermions) is factorized and thus satisfies the Yang-Baxter 
equation (with spectral parameters). 

A general consequence of integrability is the existence of a commuting family of transfer 
matrices that contains the original transfer matrix describing the discrete time evolution 
of the system. In the case of free fermions, these transfer matrices precisely "grow" Schur 
functions. Here we find in fact two families of commuting matrices [2]. It is quite satisfying 
that computing their matrix elements naturally produces one of the expressions that define 
Littlewood-Richardson coefficients, namely the coproduct formula. We thus obtain a direct 
proof of the combinatorial formula for them. 

Furthermore, the integrability strongly suggests to introduce arbitrary spectral parameters 
into the model: this corresponds to extending the original tiling model to a more general in- 
homogeneous model. This naturally produces generalizations of the Littlewood-Richardson 
coefficients which are polynomials of the inhomogeneities. We recover this way several known 
formulae [10, 16] as well as a new one for the coefficients in the expansion of the product of 
factorial (or double) Schur functions. 

The paper is organized as follows. Section 2 is a presentation of the tiling model which 
will be used throughout this paper. Section 3 presents the main ingredients in the derivation 
of the coproduct formula: the Fock space/transfer matrix formalism. Section 4 discusses 
the integrability of the tiling model and provides the proof of the main theorem used in 
section 5, which is the derivation of the coproduct formula. Finally, section 6 describes the 
inhomogeneous model and its application to product formulae for factorial Schur functions. 
The appendix briefly discusses the equivalence to a square-triangle-rhombus tiling model. 
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Figure 2. An example of tiling, using tiles a±, /3±o, 7_. 

2. The tiling model 

We provide here our own formulation of the model of tilings by decorated rhombi and 
triangles which is the basis of [10] (we actually extend it slightly by defining two additional 
tiles). 

2.1. Tiles. The model is defined by filling a domain contained inside the triangular lattice 
with tiles of the form shown on Fig. 1. More precisely, one can use either the colored lines 
inside the tiles or the symbols on the edges to decide if adjacent tiles match: that is, symbols 
must coincide on the edges, or equivalently green and red lines must propagate across edges. 
Colored lines or edge labels can be thought of as two equivalent ways to encode the four 
possible states of edges, which each have their advantages. We shall mostly use colored lines 
in what follows. The correspondence of edge labels with the notation of [10] is — <-> 0, 
+ 1, 10, 01. The shading of the tiles 7 will be explained later (cf a similar 
shading in [10]). 

Note that we have classified tiles in pairs; this is because the two tiles of type /3+, 
Po, 7_, 7 + , 70 always appear together on adjacent triangles to form rhombi, which we have 
illustrated by using dotted lines. There is some freedom however in how to reconnect the 
two tiles of type a±: an upper tile of type a_ can have on its left any number of pairs of 
tiles of type the only way this series can end is with a lower tile of type a_. Similarly 
for the tiles of type a + , with series of tiles of type /3_ on its right. In what follows it will 
be convenient to consider tilings of the whole plane. However, we shall see that with our 
"boundary conditions" (conditions at left and right infinity), all such series of tiles of type 
j3± will be necessarily finite. 

2.2. Paths. At this stage we can forget about the underlying triangles and rhombi and 
simply keep track of the paths formed by the green and red lines. Consider a horizontal line 
in the triangular lattice. Each edge can be in only three states: empty or occupied by a green 
or red line. In what follows we shall number edges using alternatingly half-odd-integers and 



4 



P. ZINN-JUSTIN 



\\ /Y VA / 

Figure 3. An example of evolution of paths from one time step to the next. 

integers to take into account the nature of the lattice. (One could of course get rid of this 
issue by applying an additional shift by a half-step say to the right, but that would break 
the left-right symmetry of the model, and we do not choose to do so here.) 

We now analyze what happens to the lines during one "time step" , that is as one moves 
(upwards) from one horizontal line to the next. Let us first ignore the tiles 7o,-,+- This cor- 
responds to the original tiling model of [19] which also occurs in [10] (in the non-equivariant 
case). Then the rules are as follows: a green (resp. red) line moves one half-step to the 
right (resp. left) if there is currently no particle of the opposite kind at this spot. The only 
situation left to consider is when lines of opposite colors are adjacent, with the green line at 
the left. Then two scenarios occur: either the green line crosses all the red lines at its right 
until it finds an empty spot, or the red line crosses all the green lines at its left until it finds 
an empty spot. This is shown on Fig. 3. 

If we add the tiles 7_ (resp. 7 + ), then green (resp. red) lines have the additional possibility 
of moving in the opposite direction as normally, on condition that the spot is free. If we add 
the tile 70, then green and red lines are allowed to simply cross each other as if they did not 
see each other. 

3. FOCK SPACES AND TRANSFER MATRICES 

We first define the notion of Fock space. The idea is that to encode the possible configura- 
tions of tilings on a given horizontal line into a Hilbert space. But first we describe another 
Fcok space which is slightly simpler (only two states per site instead of three) and will play 
an important role. 

3.1. Fermionic Fock space T. The Fock space T is an infinite dimensional Hilbert space 
with canonical orthonormal basis defined as follows. Each element |/) of the basis is indexed 
by a map / from Z + | to {— 1,1} such that there exists iV_, N + such that f(i) = —1 for 
i < AL and f(i) = +1 for i > N + . Call N + (f) (resp. iV_(/)) the smallest (resp. largest) 
such integer. 

We shall represent the — l's (resp. l's) as green (resp. red) particles or dots. 

There is a notion of "charge" which can be thought of as follows: each green particle has 
charge —1 and each red particle has charge +1. This is ill-defined because there is an infinite 
number of particles, so we need a reference state. Define \0) (the vacuum state) to be the 
state such that there are only green particles to the left of zero and only red particles to the 
right. The corresponding map from Z + \ to { — 1, 1} is the sign map. \0) has by definition 
zero charge. This way, the charge of any state |/) is given by c(f) := ^ ig z + i (/(«) — s ig n (^))- 
The charge is always an even number (we use twice the standard convention, for reasons 
that will become clear). 
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Define the shift operator S: it is defined by S |/) = |/') with f'(i + 1) = f(i) for all i. S 
decreases the charge by 2, and is an isomorphism between subspaces of different charge. 

In a subspace of given charge, basis elements can alternatively be indexed by Young 
diagrams [7] (see also [20]). The correspondence goes as follows. Rotate the Young diagram 
45 degrees, assign green dots and red dots to edges of either orientation as indicated on the 
picture: 




One can then flatten the line and obtain a configuration of green and red dots. There 
remains the arbitrariness in shifting the line, or equivalently in the charge. Here we shall 
only consider the case of zero charge, in which case the convention is that the diagonal line 
(dotted line on the picture) represents the zero. In such a way to any Young diagram A we 
associate a state simply denoted by |A), and all the basis vectors of the subspace with zero 
charge are recovered this way. In the case of the empty diagram we recover our vacuum state 
|0>. 

Finally, define for future use T+jt (resp. T-^) to be the span of the |/) such that N+(f) < k 
(resp. N_{f) > -k). 

3.2. Fock space Q of the tiling model. The Fock space Q can be similarly described as 
follows. Basis vectors |/) of Q are indexed by maps / from Z + \ to {—1, 0, +1} such that 
there exists AL, N + such that f(i) = —1 for % < iV„ and f(i) = +1 for % > N + . 

The correspondence to configurations of the tiling model described in the previous section 
is as follows: each basis vector of Q encodes a horizontal line in a configuration of tiles; thus, 
— 1 = — correspond to a green particle, to an empty spot and + 1 = + to a red particle. 
Since the model is translationally invariant the choice of an origin is irrelevant; however, 
note that successive lines have all sites shifted by half a step, which means that Q can only 
describe rows of a given parity, not both at the same time. We shall come back to this point 
below. 

Define N±(f) similarly as before: that is, iV_(/) is the location of the leftmost empty spot 
or red particle minus one half, whereas N + (f) is the location of the rightmost empty spot 
or green particle plus one half. We can also define two more numbers which will be useful: 
N-o(f) is the location of the leftmost red particle minus one half, whereas N + o(f) is the 
location of the rightmost green particle plus one half. 

There are two "quantum numbers" in Q. The first one, the charge, is defined in Q as in 
T by c(f) := X/iez+i(/(0 — sign(i)) i- e - green particles have charge —1, red particles have 
charge +1, empty spots have zero charge. The charge is an integer with arbitrary parity. 
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Figure 4. The concatenation map. 

The second quantum number, the "emptiness number" , is simply the number of zeroes: 
e(f) := #{z : f(i) = 0}. 

Intuitively, the conservation of the two quantum numbers is related to the conservation 
of the number of lines of either color in any finite region. Since we have an infinite system, 
particles can however "leak to infinity" , which results in variation of these quantum numbers. 

There is again a shift operator, denoted by S, which to |/) G Q associates |/') such that 
f'(i + 1) = f(i) for all i. S preserves the emptiness number, and decreases the charge by 2. 

3.3. From T to Q. There are two types of maps we need to define from T to Q. 

There is the obvious inclusion map from T to Q . This identifies T with the subspace of Q 
with zero emptiness number. 

The less obvious map U takes two basis elements |/_) and |/ + ) and produces |/) = 
|/_) U |/+) in Q such that 

fm = /K/-«-i) *<o 

In other words it "concatenates" the two words by discarding the right of /_ and the left of 

More generally, define for k G Z 

|/_)U fc |/ + ) = S- fc |/_)US fe |/ + ) 

This map U k is injective if one restricts to |/_) and |/ + ) such that N + (fJ) < k and N_(f + ) > 
—k, which is the only situation where we shall use U^. We thus consider from now on as 
a linear map from JF + fc ® jF_ fc to Q. 

It is an easy calculation that if |/) = |/_) U k \f + ), 

c(/) = i(c(/-) + c(/ + )) 
e{f) = \{cU-)-c{f + )) + 2k 

The image of U^, denoted by Qf rcc C Q, is exactly the span of the |/) such that N_ (f ) > 
and N +0 (f) < 0. 

Remark 1: intuitively, this second operation has the following meaning. When the sets 
of green and red particles are widely separated from each other (green ones being on the 
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Figure 5. An example of evolution of the free fermionic model. 

left and red ones on the right), then each of them behaves like a system of fermions (the 
fermionic character being the condition that there can be at most one particle per site). 

Next we shall define transfer matrices. In fact, we should say a few words on what we mean 
by a "transfer matrix" here because of the fact that we are dealing with infinite-dimensional 
spaces. A transfer matrix is defined here as a matrix, that is a collection of entries (Tf jS ) 
where / and g index the canonical basis of T or Q . It is tempting to associate to it a linear 
operator T on T or Q, such that Tf t9 = (f\T \g), but this is problematic because sometimes 
the action of T leads to an infinite linear combination of basis elements, which would require 
discussing the convergence of summations. However in all that follows, whenever we have 
two transfer matrices T and T', the product (TT') f jh = J2 g T ' f,g^'g,h om y involves finite 
sums and is therefore well-defined; so that we can safely ignore this subtlety and manipulate 
transfer matrices as operators. 

3.4. The transfer matrix of free fermions. We first define the usual dynamics for free 
fermions that leads to Schur functions, see e.g. [20]. 

The transfer matrix, denoted by Tf ree (-u), is most simply described by considering red dots 
as lines propagating (similarly as green and red lines in the tiling model). In this case the 
rule of evolution for red lines is that at each step, they can either move straight upwards 
or upwards and one step to the right on condition that no two lines touch each other. An 
example is given on Fig. 5. Furthermore, sufficiently far to the right, we impose that red 
lines go straight upwards. This way any evolution only involves finitely many moves to the 
right: we then assign a weight of u to each such move. Explicitly, (/| Tf vee (u) |/') equals the 
sum over configurations of the form of Fig. 5 where the initial configuration (at the bottom) 
is described by /' and the final configuration (at the top) is described by /, of u to the power 
the number of moves to the right. 

One remark is in order. One can of course also assign lines to the green dots and formu- 
late the rules in terms of the green lines (see [20] for details). These lines have also been 
represented on Fig. 5. The rule is that at each time step green lines can move up half-way 
then arbitrarily far to the left then up again, but in such a way that they do not touch any 
other green lines along the way. The weight of u is given to each crossing of green and red 
lines. 

Tfre^w) breaks the "particle-hole" symmetry that exchanges left and right, green and 
red lines, since the rules are clearly different for the two types of lines. One can therefore 
introduce a mirror-symmetric transfer matrix Tf TCC (u). It is defined similarly as Tf rce (w), but 
this time, the green lines are allowed to go either straight upwards or upwards and one step 
to the left. Each left move is given a weight of u. 

Both Tf rce (u) and Tf rcc (u) preserve the charge. 

Finally, we have the following important formulae: 
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Lemma 1. 

n 

SX/n(Ul, ■ ■ ■ , U n ) = (fl\ Jl T frce (Wi) |A) 

i=l 

is the skew Schur function associated with the Young diagram X. For \x = 0, 

n 

s x {ui, . . . , u n ) = (0| JJ T bee (ui) | A) 

i=l 

is the Schur function associated with the Young diagram A. Similarly, 

n 

Sxt/^t^Ui, ...,Un) = (fi\ Y[T lrcc (ui) | A) 

1=1 

where X T is the transpose of X, and in particular 

n 

S x t(u U . . . , Un) = (0| Yl Tfrcc(^i) |A) 

i=l 

The most general formula is 

m n 

s\/n( u i> • • • j u m /vi, . . . , v n ) = (n\ Y\ Tf ree (wi) T bee (vi) |A) 

i=l t=l 

t/iat t/ie super symmetric skew Schur function, which leads for \x = to the usual super- 
symmetric Schur function 

m n 

S X (U U . . .,U m /vx, ...,V n ) = (0| JjT frcc (Mi) JjTfrcc^) I A) 

1=1 1=1 

Remark 2: as will be apparent in the proof, the expressions in the lemma are independent 
of the ordering of the products, which is why we left them unspecified. This implies the 
commutation relations 

[T frcc (n), T free (»] = [T frcc (n), T frcc (t;)] = [T hee (u),T hee (v)] = 

which are also a consequence of the more general results of the next section. 

Proof. There are several simple proofs of this standard result. Note first that taking products 
of n transfer matrices amounts to stacking together n rows made of the paths defined above. 
One proof involves a bijection between these paths and the appropriate tableaux that one 
uses to define supersymmetric Schur functions. Another proof, which we sketch here, is to 
use the Lindstrom-Gessel-Viennot (LGV) formula [14, 6]. We apply it to the green lines to 
the right of iV_(A) (those to the left necessarily go straight). There are exactly X[ such lines, 
where X[ is the number of non-zero rows of A. This leads us to compute the evolution for a 
single line from initial location k to final location k + j. Noting that the rules of evolution 
are translationally invariant, one can introduce a generating function h^(x) = ^j^hjxi 

for the time step corresponding to Tf ree (itj) and e^(x) = Y^j>o e j x '' f° r ^ ne time step 
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corresponding to Tf ree (i>j). hj and ef 1 are the numbers of ways to move j steps to the left 
for a single green line and a single time step, so we immediately find 

h {i \x) = (1 - xui)~ l e {i) (x) = 1 + xvi 

Composing the transfer matrices amounts to multiplying the generating series, so we find 
the evolution for a single green line to be given by the generating series 

h( i \ - nr=i( i +xvj) 

which is exactly the generating series of the supersymmetric analogues of completely sym- 
metric functions i.e. Schur functions hj corresponding to one row: 

h(x;ui, . . .,u m /vi, ...,v n ) = y^hj(ui, . . .,u m /vx, . . .,v n )x J 

Applying the LGV formula produces the Jacobi-Trudi identity for supersymmetric Schur 
functions 

sa/ m (ui, . . .,u m /vx, . . . ,v n ) = det (hx.- IH+ i-j(ui, . . .,u m /vx, . . ..Vn))^. .^ 

□ 



3.5. The transfer matrix of the tiling model. Before defining the transfer matrix of the 
tiling model, one must discuss the problem of the numbering of successive rows. Since the 
edges are shifted by a half-step during one time unit, there are two symmetric choices: shift 
everything one half-step to either the right or the left (either before or after the evolution, 
since it is translationally invariant). The resulting transfer matrices are called respectively 
T + and T_. (/| T± \g) is then defined as the number of possible ways paths can move from 
the initial configuration encoded by g to the final configuration encoded by / in one time 
step according to the rules of evolution described in section 2.2, using only the 2x5 tiles a 
and P of Fig. 1. In principle one could assign weights to the different tiles but we shall not 
need to do so here. Note the relation T + = ST_. In the case of a two- row evolution, one 
can introduce T 2 := T + T_ = T_T + , which has the advantage that it restores the left-right 
symmetry. 1 

Note that sufficiently far to the left, there are only green lines and these necessarily move 
one half-step to the left. On the contrary, far to the right, one has red lines that move one 
half-step to the right. This observation allows us to conclude that T± changes the charge by 
=1=1 and increases the emptiness number by 1, or equivalently T 2 preserves the charge, but 
increases the emptiness number by 2. 

We now prove some additional properties of T 2 . 
Lemma 2. For any pair of basis states f and g, 

(/| T 2 \g) ± M (/) > N- (g) + 1 and N +0 (f) < N +0 (g) ~ 1 

1 Despite the notation, T 2 is not the square of an operator on Q; if one insisted that it be so, one would 
have to define T as acting on Q ®Q, with T = ( ^ T ^ J . 



10 



P. ZINN-JUSTIN 



Proof. Since particles of the same color never cross, the leftmost or rightmost particles remain 
the same during time evolution. The lemma then follows from the fact that red (resp. green) 
particles move at least one half-step to the right (resp. left) at each time step. □ 

Lemma 3. If\f) = \f-)U\f+) G G bee , with JV+(/_) < and N-(f+) > 0, then T 2 |/) e S free; 
and 

T 2k \f) = \f_)u k \f + ) Vkez + 

Proof. Such a state |/) has the properties that all green lines are to the left of red lines, so 
no crossings ever occur. Therefore, the green (resp. red) lines move one half-step to the left 
(resp. right) at each time step. This is all that the lemma says. □ 

Lemma 4. Call p = m&x(N +Q (f), — N_ (f)). Then T 2p \f) G Gf Tee and there exist unique 
coefficients c gh such that 

T 2k \f) = Y, c l,h\9)Uk\h) Vk>p 

9,h 

Intuitively, this lemma says that no matter what the initial state is, eventually all possible 
crossings will take place and we shall be left with a linear combination of states which are 
all such that all green lines are the to the left of red lines. 

Proof. This is combination of the two preceding lemmas. Let |/) be a basis state, and set 
p = max(iV + o(/), — iV_o(/)). Then according to lemma 2, T 2p \ f) is a linear combination of 
basis states \g) such that N-o(g) > and N +Q (g) < 0, i.e. T 2p |/) e feee- 

By definition of Ghee, this implies that there exist uniquely defined coefficients c gh such 
that T 2p \ f) = J2 g h c gh \d) 4> 1^)' where the summation is restricted to couples (g, h) such 
that N + (g) < p, N_(h) > —p. Set all other entries c^ h to zero. The formula of the lemma 
then follows by application of lemma 3. □ 

Corollary. If |A) G T C Q is a state with zero charge and zero emptiness number and 
p = max(iV + (/), —N_(f)), then there exist unique coefficients c^ u such that 

T 2fc |A) = ^<J/i)U fe |z/) \/k>p 
where the sum is over pairs of Young diagrams. 

Proof. Such a state |A) has zero charge and zero emptiness number. The formula of lemma 
3 implies that all pairs (g, h) that contribute to the sum satisfy = c(A) = \{c{g) + c(h)) 
and = e(A) = \(c(g) —c(h)), so that they have zero charge themselves. They can therefore 
be indexed by Young diagrams. □ 

3.6. The two families of commuting transfer matrices. We define here two one- 
parameter families of transfer matrices T±(u), which are closely connected to T±. 

Naively, one would like to define T± as transfer matrices which, in addition to the usual 
tiles a and (5, allow one more type of tiles ■y±, but at a cost of a certain weight x. However, 
this turns out to be inconvenient for the boundary conditions we have in mind, so we use 
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the following definition instead, using the inverse u = x~ l of the natural parameter x (which 
will reoccur in section 4). T±(u) is defined similarly as T±, but with three new ingredients: 
(i) we allow the extra tiles of type 7± and (ii) we give a weight of u to each pair of tiles of 
type a±, and to each pair of tiles of type f3± (all the signs being the same as that of the 
transfer matrix), (iii) we impose the following condition at infinity: sufficiently far to the 
left and to the right, the effect of T_(iz) (resp. T + (u)) must be to push either type of lines 
one half-step to the right (resp. left). Taking into account that this half-step is absorbed into 
the definition of the matrix, we find that at infinity T±(u) behaves like the identity. This 
ensures that only a finite number of tiles of type a± or f}± ever occur in the evolution w.r.t. 
T±(tt). Thus each entry of T±(u) is a polynomial in u. 

The boundary conditions also imply that T±(u) preserves the charge and the emptiness 
number. 

We now list some properties of T±(u). 
Lemma 5. For any pair of basis states f and g, 

(f\ T+(u) \g)^Q M (/) > N- (g) - 1 and N +0 (f) < N +0 (g) 
(f\ f _(«) \g)^0 iV_ (/) > N_ (g) and N +0 (f) < N +0 (g) + 1 

Proof. Same proof as lemma 2, but this time green particles can move one half-step to the 
right for T_ and red particles can move one half-step to the left for T+. □ 

Lemma 6. |/) = |/_) U |/ + ) G Q iicc , then 

T_(«)|/) = T faB («)|/_>U|/ + > 
f+(ti)|/> = |/_>UT faB («)|/ + > 

Proof. By inspection. If there are only green particles to the left of and red particles to the 
right of 0, no crossings can occur. For T_(it), this leaves two possibilities for green particles: 
going left one half-step (pair of tiles a J) with a weight of u or going right one half-step (pair 
of tiles 7_) with a weight of 1; and only one possibility for red particles, going right one-half 
step (pair of tiles a + ) with a weight of 1. Adjusting the locations by shifting everything one 
half-step to the left, we obtain exactly Tf ree (ii) for green particles and no evolution for red 
particles. The reasoning is the same for T + (u). □ 

Note the similarity with lemma 3. However, an important difference with the situation of 
lemma 3 is that it is not true that |/) G Gf Tce implies T±(u) |/) G <?f rcc , so in general one 
cannot iterate the argument. 

Remark 3 (followup of Remark 1) : intuitively, this means that when the sets of green and 
red particles are widely separated from each other then our transfer matrices make them 
evolve like free fermions. When green and red particles mix, it is not the case any more: in 
this sense what we have is a system of two species of fermions interacting with each other. 
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Lemma 7. The transfer matrices T±(u) leave T C Q stable, and: 

\f) = T free (w) |/) V|/)G^ 

f +(U) |/) = Tfree(u) |/) V|/)6^ 

Proof. T viewed as a subspace of Q is simply the subspace of zero emptiness number, and 
the T±{u) preserve this number, so they leave T stable. 

The rest of the reasoning is again by inspection. Consider the action of T + (u) on a state 
with no empty spots. Remembering that the final state should also have no empty spots 
(this is in fact a consequence of the boundary conditions, as already explained), we conclude 
that tiles of type a + are forbidden (consider the leftmost tile of type upper a + ; the only 
allowed tile left of it is the empty tile fi ). Similarly tiles of type /3_ are forbidden (sequences 
of such tiles always end with tiles of type a+). All that we are left with is half-steps to the 
left (tiles a_ and 7+) and crossings of the type /?+, that is one green line crossing a series 
of red lines. But up to an overall half-step to the left, these crossings are exactly those that 
occur between red and green lines in the free fermionic model, compare Figs. 3 and 5. The 
weight of u is given to each each pair of j3 + tiles, that is to each crossing, which is the same 
weight that is given in the free fermionic model. 

A similar reasoning works for T_(it) and Tf ree (w). □ 
Theorem 1. All transfer matrices T_(u), T + (v), T commute: 

[f e (u), f e ,(v)} = [f e (u), T e ,] = e,e'e {+, -} 

This is the central result of this section, embodying the integrability of the model. Its 
proof is the subject of next section. 

4. Yang-Baxter equation and proof of the commutation theorem 

We now discuss the integrability of the tilings model, that is the various forms of the 
Yang-Baxter and how they imply the commutation relations of theorem 1. 

4.1. -R-matrix and Yang— Baxter equation. The first idea is to decompose the action of 
the transfer matrices as products of elementary blocks, called i?-matrices. It is convenient 
in this section to go back to the labels +,—,0,0 on the edges because using it, the tiles 
are explicitly Z 3 -invariant. We can then choose as elementary blocks two adjacent triangles. 
They have three possible orientations, which correspond to three possible rhombi. 

Let us choose for example the orientation as shown on Fig. 6, and allow the tiles a±, /3±o, 
7 + , with the weights indicated below, x is the spectral parameter. In fact the only rule is that 
pairs of shaded tiles get a weight of x. The weights look different but are actually equivalent 
to the weights of T + (v), with v = 1/x (this will be explained in the next section). The 
colored lines are only decorative here and the focus is on the edge numbers. We can encode 
this information into a3x3x3x3 tensor as follows. R^fc^x), where each index lives 
in { — , +, 0}, is the weight of the rhombus with edges k, I read counterclockwise starting 
from the left edge (or the right edge, the tiles being 180 degrees symmetric). 
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FIGURE 6. The i?-matrix in one of its three orientations. 
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R(x) is usally considered as a matrix by associating to it the map from C 3 ® C 3 (corre- 
sponding to the left and bottom edges) to itself (corresponding to the right and top edges); 
explicitly, 



(4.1) 
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This point of view is not particularly useful for our purposes because it breaks the Z 3 
rotational symmetry by forcing to distinguish "incoming" and "outgoing" edges. 

Rotate now every tile 120 degrees clockwise: we obtain the tiles a±, (3±o, 7_ (which are 
exactly the tiles of T_(v)). This time it is 7_ pairs (again, the shaded tiles) which get 
a weight of x. Finally, there exists a third possible orientation of rhombi, obtained from 
the first set of tiles by 120 degrees counterclockwise rotation - and correspondingly, a third 
shaded tile which has not been used so far, namely 70. We shall use the following graphical 
notation for these three types of /^-matrices: simply depict them with rhombi (with thick 
edges) and put the spectral parameter x inside. The convention is that when several rhombi 
are pasted together, the internal edges are free i.e. summed over while the external edges are 
fixed. 

We now have the key proposition, which is the Yang-Baxter equation: 

Proposition 1. Let three variables x, y, z satisfy x + y + z = 0. Then, the following equality 
holds 





for any values of the external edges while summing over the matching tiles inside. 

Proof. The proof is just an explicit computation. We provide a particularly compact version 
of it, which was suggested to the author by A. Knutson. 

Note first that if the equation is true for a given sequence of external edges, then it is true 
for any cyclic rotation of it: indeed a rotation of 120 degrees amounts to cyclic permutation 
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of the variables, and a rotation of 180 degrees amounts to exchanging l.h.s. and r.h.s. (the 
tile weights being invariant by 180 degrees rotation). 

Next, observe that if no shaded tiles appear in the equality, then it is trivial - the different 
types of rhombi are only distinguished by the additional shaded tiles which carry a non-trivial 
weight. If there are such pieces, then necessarily one of the external sides of the tiles must 
be a sequence of — , + (read counterclockwise). This leaves only three possible sequences 
of external edges, up to cyclic rotation: (i) —,+,—,—,+,—, (ii) —,+,+,—,+,+ and (hi) 
— ,+,—,+,—,+. The first two sequences are invariant by 180 degrees, which allows to prove 
the equality without calculation, since as already mentioned the 180 degrees rotation of the 
l.h.s. is the r.h.s. The last case is the only interesting one: we find the identity 




(or its 180 degrees rotation), that is x + y + z = 0. □ 



4.2. The RTT relations. We use the notation T± 6 (x) for the transfer matrix that creates 
one extra row made of the usual tiles a, (3 and 7±, followed as usual by a ±|-step, in such 
a way that sufficiently far to the left the left sides of the rhombi are all a and sufficiently 
far to the right the right sides of the rhombi are all b. The weights are those of R, divided 
by a factor hi± if say the value of the upper edge is equel to ±, which makes sure that 
the tiles have weight 1 sufficiently far to the left and right. Explicitly, K a + = R aj+ a + and 
K°2 b = R_ j j, : _ j 6. In practice we shall only consider here the cases T±' + , T^'~ and Tt' + , in 
which case we recover. 

T ± = T±' + (0) T+(u) = T;'-(0 T_(u) = Tl'+(0 

We now apply the standard argument of repeated application of the Yang-Baxter equa- 
tion ("unzipping") to write the RTT relations, which allow to show various commutation 
relations. 

Let us fix two basis elements |/) and \g) in Q and consider matrix elements of the form 
(f\T°^(x)T'z d (y) \g). The important observation is that to the left of min(iV_(/), iV_(g)), 
one has series of identical tiles. Thus no information is lost by truncating the state anywhere 
left of min(iV_(/), N_(g)), and the product of weights of the discarded left tiles is 1. Similarly, 
to the right of ma,x(N + (f) , N + (g)), everything is frozen and we only have weights of 1. 

Note that we already know this to be the case "sufficiently far to the left", but this is 
not good enough for our purposes: we need a bound that is uniform w.r.t. the intermediate 
state. For example if a = b = — and c = d = +, one can check that the infinite sequence 





on the left side and of 
lower edges change their value. 
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In the case a = b = — and c = —, d = +, there is no change to the right of max(iV + (/), N + (g)), 



+ 



where we find again sequences of 




. However, to the left of min(iV_(/), N_(g)), we have 



+ 



a new situation where green lines move one full step to the left, that is sequences of 



The two other cases of interest to us can be treated similarly. 

Thus, if we pick any region containing the interval [min(iV_ (/), N_(g)), max(iV + (/), N + (g))], 
we can write that 

/ 

Wn^r (/i nvjr/fo) \ 9) = ) ,),),) -) ,),),/ 

9 

where g (resp. r) denotes the number of green (resp. red) lines in the finite domain, which is 
well-defined if b = — (resp. c = +) because no green (resp. red) line can cross the boundary 
with such boundary conditions (and these are the only cases where the factors k are not 1). 

At this stage, we can insert the i?-matrix at the left, with spectral parameter z = —x — y: 

/ 

E^.,( : )K'na' (f\n\xm*(y) | 9 > = {) ) ) ) ■) ) ) ) 

9 

We apply repeatedly proposition 1 all the way to 



/ 



J2^c,a',cA^tn^ d ) 9 (f\ Tf (x)T^(y) \g) = 

a,c 




9 

= J2 R «^»(- x - y)(4' b 'n^ d ') 9 (f\ Tt' d ' ( y )T a ;> b ' & \ g ) 

b',d' 

By varying a', b, c', d one gets in principle 81 identities relating transfer matrices with 
various boundary conditions. We are only interested in the four cases a' = — , b G {+, — }, 
c'g{+,— }, d = + which lead to the following 

Theorem 2. We have the commutation relations 

[Ty-(x),Tj + (y)] = [T^(x) 1 Tt' + (y)} = [T~' + (x) , TZ' + (y)] = [T+> + (x), T + _' + (y)} = 

This theorem implies that [T+(«), T-(u)] = 0, as well as [f e («), T e /] = 0, e, e' G {-, +}. 
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Figure 7. Graphical representation of the Yang-Baxter equation with one 
matrix R*(x) and two matrices R(y), R(z). Here x = z — y. 

4.3. Other cases. In order to complete the proof of theorem 1, we must show the com- 
mutation [T±(w), T±(u)] = (same sign). This requires another kind of i?-matrix, as the 
geometry suggests, cf Fig. 7. We do not provide the details since these commutation relations 
are not needed anywhere. Let us simply give the expression of the new R- matrix: 
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1 / 



The rest of the proof is identical (unzipping argument) and we obtain the following 
Theorem 3. We have the commutation relations 

[T^(x),T+>-(y)\ = [Tt' + (x),Tt' + (y)\ = [T7j+(x), T~> + (y)] = [Tj + (x),Tj + (y)} = 

Note that this theorem does not say that T+'~ and T+' + commute. 

5. Littlewood— Richardson coefficients from the coproduct 
5.1. The coproduct formula. 

Theorem 4. The coefficients c^ v in the corollary of lemma 4 ore Littlewood-Richardson 
coefficients. 

Proof. Start from the supersymmetric Schur function (lemma 1) 

m n 

S A (ui, . . .,U rn /v 1 , ...,V n ) = (0| JjTfree^) JjTfree^i) |A) 

i=X i=l 

Considering |0) and |A) as states in Q, we immediately find by applying lemma 7: 

m n 

s x (u 1 ,...,u m /v 1 ,...,v n ) = (0| JJf + (M<) JjT_(ui)|A) 

1=1 1=1 



LITTLEWOOD-RICHARDSON COEFFICIENTS AND INTEGRABLE TILINGS 17 




AT- (A) N+(\) 



Figure 8. An example of configuration and the triangle in which all the 
crossings take place. 

Given a non-negative integer k, consider the evolution backwards in time (i.e. going down- 
wards instead of upwards) starting from (0| l_4 (0|. It is elementary to check that after 2k 
steps, one gets (0|, or in other words (0\ = (0| U k (0| T 2fc . This results in 

m n 

s a (mi, . . . , u m /vi, ...,«„) = (0| U fc (0| T 2fe JJ f n f _(«() I A) 

i=l i=l 

We now use the commutation of the transfer matrices (theorem 1) to rewrite it 

m n 

m, . . . , o = <0i u fe (0| n t n f_( Vl )T 2k i a) 

i=l t=l 

Assume fc > max(A r + (A), — AL(A)), and apply the corollary of lemma 4: 

s A (ui,...,u m /ui,...,^) = ^c^(0| U fc (0| JJf+(Mi) JjT-(ui) U fc |z/) 

/Lt,i' i=l i=l 

We now wish to apply lemma 6. \fj,) \v) G feee, but as already mentioned T + (u) may 
push red particles to the left and T_(u) may push green particles to the right. Let us thus 
choose k > max(iV_|_(A), — A r „(A)) + n/2; according to lemma 5, this ensures that lemma 6 
can be applied repeatedly. Separating into two bra-kets, we get: 

m n 

sx(ut, . . .,u m /vi, ...,v n ) = y^c^ v (0\ T frco (u;) \v) (0\ T bee (v i) \fi) 

H,u i=l i=l 

Finally, apply the basic free fermionic identity (lemma 1): we find 

Sx(u u . . .,U m /v U ...,V n )= 7~]c^ S„(U U . . .,U m )s^r(v U . . . , V n ) 

Remembering that the l.h.s. is the supersymmetric Schur function associated to A, we con- 
clude that this equality defines uniquely the as Littlewood-Richardson coefficients. □ 
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5.2. Back to the triangle. All that has been described so far uses an infinite-dimensional 
Fock space of configurations. This helps in formulating the results more elegantly. However, 
for any given Young diagram A, all the evolution of the state |A) with respect to T takes 
place in a finite part of the plane; more precisely, a triangle, which leads us back to the more 
standard formulation, as in [10], in terms of puzzles. 

First, let us define how to read Young diagrams from the boundaries of a finite domain 
(this is a finitized version of the transformation of section 3.1). Let n and k be two non- 
negative integers such that k < n. Consider a sequence of n successive edges which can be 
either — or +, such that there are k "— " and n — k "+". Fix an orientation of the edge: 
then reading the values of the edges following the orientation, one obtains a binary string 
of length n. To it one associates bijectively a Young diagram contained inside the Young 
diagram k x (n — k), that is the rectangle of height k and width n — k, as follows: each 
— corresponds to a step up, and each + to a step to the right, starting from the lower left 
corner of the rectangle. For example, 



Note that depending on the direction of the boundary, a — will be drawn as a green 
particle or an empty spot, and a + will be drawn as a red particle or an empty spot. 

Let us now fix an initial state given by a Young diagram A and consider the possible 
evolution of the system. The situation is illustrated on Fig. 8. The smallest possible triangle 
in which all crossings between green and red particles take place has as its lower side the 
interval between N_(X) and N + (X). So we can set n = iV + (A) — iV_(A), in which case k is 
the height of A (number of non-zero parts). Of course one could make the triangle bigger, 
which illustrates the stability property of puzzles. Now the configuration outside the triangle 
is uniquely fixed by the locations of the lines at the two upper sides of the triangle: green 
(resp. red) particles move uniformly to the left (resp. right). In other words, no information 
is lost by restricting to the triangle and conversely, any configuration inside the triangle can 
be extended to the outside in a unique way. The "asymptotic" states \i and v which describe 
the sequences of green particles and empty spots to the left and of red particles and empty 
spots to the right can also be read off the two upper sides of the triangle in the way described 
in the previous paragraph, the binary strings being read from left to right. Combining these 
observations, we conclude that v also counts the number of fillings of the triangle with 
fixed boundaries, where A is encoded by the bottom edge, \i by the left edge and v by the 
right edge, always reading from left to right. These fillings are called 'puzzles. For example, 
on Fig. 8, one reads 



+-++-+ 




A 



, n = v = 



6. Equivariance or the introduction of inhomogeneities 



It is natural to try to generalize the construction of section 5 by putting arbitrary horizontal 
and vertical spectral parameters (inhomogeneities) at each site. As already observed in 
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[4, 12, 3] in a different but similar setting, this amounts geometrically to going from ordinary 
cohomology to equivariant cohomology - in the present case, of the GraBmannian. 

There are however some complications. The first one is that the infinite system viewpoint 
used so far is not particularly convenient to handle the functions that appear in the inho- 
mogeneous situation, that is factorial (or double) Schur functions. We shall therefore work 
in a finite domain instead. 

The second one is that the coproduct formula does not obviously generalize to the inho- 
mogeneous case. In fact, a coproduct formula for double Schur functions has only recently 
been discovered [17] and the description in terms of objects analogous to puzzles is unknown. 
We shall prove instead product formulae. Note that even specialized to the non-equivariant 
case, the formulae of this section will thus be distinct from those of section 5. 

The basis setup is to consider that there are (oriented) lines of spectral parameters propa- 
gating on the dual lattice of the lattice of tiles of i?-matrices, in such a way that the spectral 
parameter attached to the tile is the difference of the two spectral parameters crossing it: 

y 

4 

L 

In all that follow we fix integers n and k, k < n, and use the same correspondence between 
binaray strings and Young diagrams inside the rectangle k x (n — k) described in section 5.2. 

We now define the two basic "building blocks" with which we shall produce non-trivial 
identities. 

6.1. Factorial Schur functions. What we really want to define is the double Schubert 
polynomial of a GraBmannian permutation. It will however coincide with the notion of 
factorial Schur function. 

Given a Young diagram A inside the rectangle k x (n — k), let us define the factorial 
Schur function S\(xx, . . . , Xk\yi, • • • , y n ) graphically as on Fig. 9, that is as a tiling of a n x k 
rhombus with prescribed boundaries: A is encoded into the edges of the left side, read from 
bottom to top, the bottom side is full of green particles and the other two sides are empty 
(in all the figures of this section, we omit entirely drawing empty spots, and use half-colored 
circles to indicate that the spot can be either empty or occupied by a particle of the given 
color, depending on the Young diagram it encodes). The spectral parameters are the Xi and 
the i/i, from left to right and bottom to top. 

An important property is the following: 
Proposition 2. s\(xi, . . . , Xk\yi, ■ ■ ■ , y n ) is a symmetric function of the Xi, i — 1, . . . , k. 

Proof. Fix % G {1, . . . ,n — 1} and apply the matrix R*(xj + i — Xi) from section 4.3 to the 
bottom edges i and z'+l. Noting that R* 3 - = 5i-5j- and R*,- o,o = di,o$j,o, we can use the 
usual unzipping argument (repeated application of the Yang-Baxter relation represented on 
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Figure 9. Our definition of factorial Schur functions, and an example in 
various graphical representations. The numbering of the Xi is reversed in the 
tableau representation. 

Fig. 7) to move the matrix R*(xj +1 — Xi) all the way to the top and then remove it. The 
result is the same picture as we started from, but with x^ and x i+ i exchanged. □ 

On Fig. 9, an example of configuration is provided, as well as two alternative representa- 
tions of it. The first one is the pipe dream of the Grafimannian permutation associated to 
A [5, 9]. Our picture can be deformed into the pipe dream picture by rotating 60 degrees 
clockwise and then distorting it slightly to make the 60 degrees angle a right angle. This way, 
the green lines are exactly the trajectories of the lines above the descent of the permutation. 
The other lines can be recovered unambiguously. The weights are now given to the crossings 
(between a red line and a green line) and take the form Xi — yj, where i is the row number 
and j the column number, both counted starting from the corner. 

The final representation is a semi-standard Young tableau of A, with the alphabet {1, . . . , k}. 
Starting from the pipe dream picture, each column of the tableau corresponds to a red line 
counted from left to right, and each number in the boxes of the column indicates the row 
at which the red line crosses a green line. An important difference is that the rows are 
numbered from bottom to top; but since the factorial Schur function is a symmetric function 
of the Xi, this is irrelevant, and we may as well assign to the tableau a weight equal to the 
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Figure 10. Definition of an MS-puzzle (k = 2, n = 4). 

product over boxes of xt — Ut+c, where T is the number in the box (i.e. the row number in 
the pipe dream picture) and c is the content of the box (column minus row of the box), in 
such a way that T + c is precisely the column number in the pipe dream picture. 

6.2. MS-puzzles and equivariant puzzles. We first define an MS-puzzle (this terminol- 
ogy is borrowed from [10]). It is simply given by Fig. 10. Three sides encode Young diagrams, 
whereas the fourth side, the upper right one, is simply a sequence of + n ~ k — fc , all read from 
bottom to top. We denote this object by e v x „(yi, . . . ,y n ; z±, . . . , z n ). 

We shall need the following series of lemmas: 

Lemma 8. In the top half of a MS-puzzle, green lines always go in straight lines (up and to 
the left). 

Proof. By inspection. Since the red lines are stacked at their rightmost positions, there is 
no available free spot for a green particle to cross the red particles by moving straight to the 
right. □ 

Lemma 9. Assume that z n+ i_i = yi, i = 1, . . . , n. Then the top half of the MS-puzzle is 
'frozen" i.e. it has a unique configuration (which has weight 1); the edges on the horizontal 
diagonal of the MS-puzzle reproduce, read from left to right, the Young diagram ji. 

Proof. If z„+i_i = yi for all i, then the rhombi sitting on the horizontal diagonal of the MS- 
puzzle have a zero spectral parameter. This implies that the shaded tiles 70 is forbidden, or 
equivalently that the edges crossing horizontally these rhombi cannot be (red and green 
particles on top of each other). Since there are k green and n — k red lines incoming, there 
cannot be an empty spot either (edge 0). So these edges can only be + or — ; according to 
lemma 8, the green particles on it are at the same locations as on the upper left edge, and 
therefore the red lines are also fixed (they must move one step to the left each time they 
cross a green line) and occupy the complementary set on the diagonal. In other words, the 
diagram \i is reproduced on the diagonal. No shaded tiles are used in the top half, so its 
weight is one. □ 
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Figure 11. From an MS-puzzle to an equivariant puzzle. 

With the hypothesis of the last lemma, the top half of the MS-puzzle can be removed 
since it is fixed and has a weight of 1. The result, after 180 degree rotation, is called an 
equivariant puzzle, cf Fig. 11. 

If we call A the weight of an equivariant puzzle, that is the tiling of a triangle with sides 
v (bottom), /i (left), A (right), all read from left to right, with the same tiles and spectral 
parameters as an upside-down MS-puzzle with z n+ i_i = Hi, we have 

e \,fj,( Z ni ■ ■ ■ i z \ \ z \i ■ ■ ■ i z n) = c p,a(^1' • • • ' Z n) 

or equivalently c^ x (zi, . . . , z n ) = e^ p (z n , . . . , z±; z±, . . . , z n ), where A is the complement of 
the Young diagram A inside the rectangle k x (n — k) after a 180 degree rotation, which 
corresponds to reading the binary string of A from right to left. 

Lemma 10. There are no equivariant puzzles contributing to c u or if fi ^ v, and one 
(which has no shaded tiles) if /i = v, so that: 

C 0,^t( 2 'l' • • • ' Z n) — Cfi^i.Zli ■ ■ ■ i Z n) — 

Proof. Considering the equivariant puzzle as the bottom piece of an MS-puzzle, apply lemma 
9 to either the mirror image w.r.t. the horizontal axis of the MS-puzzle with colors inter- 
changed, or its 180 degree rotation. □ 

Finally, note that if all spectral parameters are equal to zero, we find by definition 
c a,m(0' • • • ; 0) = c a,/J which implies e^ M (0, . . . , 0; 0, . . . , 0) = Cp A . Rotating a (non-equivariant) 
puzzle 60 degrees counterclockwise results in another puzzle, so that A = c\ . To summa- 
rize, 

c^(0,...,0) = e^(0,...,0;0,...,0)=c^ 

6.3. The Molev— Sagan problem. The Molev-Sagan (MS) problem consists in expanding 
the product of two factorial Schur functions with the same first set of variables as a sum of 
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Figure 12. Graphical solution of the MS problem (n = 5, k — 3). 



factorial Schur functions. It was first solved in [16] in terms of barred tableaux and then 
reformulated in [10] in terms of MS-puzzles. We now rederive it in our framework. 

Let us consider the formal equality of Fig. 12. On each side of the equality, we have 
configurations where all the external edges are fixed. The order of the parameters x/i and 
Zi is important; note in particular the dashed line which corresponds to the difference of 
spectral parameters vanishing. 

On the upper left side, we have a Young diagram /i encoded by a binary string of green 
particles and empty spots in the region of parameters yi, and emptiness above in the region 
of the Zi. On the lower left side, we have a Young diagram A encoded by a binary string 
of green particles and empty spots. Both diagrams are read from bottom to top. The right 
sides have k green particles at their highest possible location and n — k red particles at their 
lowest possible location. The top and bottom sides are full of green particles. Furthermore, 
we assume that the sum of widths ji\ + Ai of the Young diagrams /i and A does not exceed 
n — k. This can always be achieved by choosing n sufficiently large. The justification of this 
assumption will become clear below. 

In order to go from one side to the other side of the figure, one simply applies repeatedly 
proposition 1 (Yang-Baxter) starting from the center and "piling up boxes". Since our 
spectral parameters are differences of the parameters of crossing lines, the sum of spectral 
parameters around a hexagon is indeed zero. 
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Figure 13. An example of tilings contributing to each side of the equality of Fig. 12. 

We shall now examine what the consequences of this equality are. It is simpler to work 
on an example, as on Fig. 13. For ease of discussion, we have labelled the various regions 
inside the hexagon and marked in blue the "interesting" ones, that is those in which there 
is some freedom left - the other regions are entirely "frozen" once the boundary conditions 
are fixed. 

Let us start with the left hand side. In region A, green particles starting from the bottom 
can only go straight to the left, so that red particles must go horizontally. In region B, red 
particles move freely to the right. In region C, we note that the spectral parameters are of 
the form of the hypothesis of lemma 9, so that we can apply lemmas 9 and 10 to conclude 
that the whole region is frozen. The result is that the upper left side of region C reproduces 
the Young diagram A read from top to bottom. We now recognize region D as the 180 degree 
rotated picture of a factorial Schur function (cf section 6.1). Since all tiles and weights are 
180 degree rotation invariant, we find that this region contributes s\(x\, . . . ,Xk\zi, . . . , z n ), 
paying attention to the order of the parameters Zi (we have however reordered the xi because 
of proposition 2). Region E is the picture of a factorial Schur function with the standard 
orientation, so it contributes s^xi, . . . ,X).\yi, ■ ■ ■ ,y n ). Overall, the left hand side equals 
s\(xi, ■ ■ -,Xk\zi, . . . ,z n )s li (x 1 , . . .,x k \yi, . . .,y n ). 

Now let us compute the right hand side. In regions G, H and I, it is easy to check that all 
particles can only move straight in the way indicated on the figure. Region F needs to be 
dealt with carefully. We want to make sure that the green lines starting from the bottom exit 
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region F through the left side. This requires to look ahead into the region J. Let us consider 
the rightmost green particle. If we number the upper left side of region J from bottom to 
top, then its endpoint is precisely k + fa (the length of the first row of jj). This means that 
the location (counted from bottom to top) of the rightmost green particle when it crosses the 
diagonal line starting at the bottom of the junction of J and F cannot exceed k + fa plus the 
number of steps to the left it can make in region J. We claim that this number is Ai: indeed 
among the n — k red particles starting from the lower left side of J, only Ai are allowed to 
cross without making a step to the right. The last n — k — Ai, which are the topmost ones, 
have to make a step to the right each time they cross a green particle in order to reach their 
final destination to the right of the upper right side of F. The result is that this location is 
at most k + fa + Xi, which by assumption is less or equal to n. Thus the rightmost green 
particle, and therefore all others, exit through the left side. We now recognize region F as 
the factorial Schur function s„(xi, . . . , XkllJij ■ ■ ■ lUn), where v is the Young diagram, read 
from bottom to top, that encodes the green lines at the boundary between F and J. And J 
is precisely a MS-puzzle (cf section 6.2), which contributes e\ Jyi, . . . ,y n ; Zi, . . . , z n ). So the 
right hand side equals J2u e x,^^ ■ ■ • > Vn, %i, ■ ■ -,z n )s u (xi, . . .,x k \yi, . . .,y n )- 
Finally, we find the desired equality 

s A (xi, . . . ,x k \z!, . . . ,z n )s, J ,(x 1 , . . . ,x k \yi, ...,y n ) 



The summation over v is only on Young diagrams inside the rectangle k x (n — k), which is 
related to the fact that (i) we have only k variables x\ and (ii) the sum of widths of [i and 
A satisfies fa + Ai < n — k. If one sets all yi and to zero, then the equality becomes 



which is the product formula characterizing the usual Littlewood-Richardson coefficients. 

6.4. Alternate solution of the Molev— Sagan problem. Interestingly, there is a small 
variation of the construction of the previous section, which produces another, not obviously 
equivalent, solution of the Molev-Sagan problem. Since it is based on the same principle, 
we shall only sketch the derivation, based on Figs. 14 and 15. 

By inspection, one easily finds that regions A, B, G, H, I are frozen, with lines going 
straight as indicated on the figure. Region C is strictly speaking not frozen; in fact, we 
recognize an upside-down (180 degree rotated) MS-puzzle, so that its weight is equal to 
e\ {z\, . . . , z n ; Z\, . . . , z n ) where A' is the upper left side of C read top to bottom. This, ac- 
cording to the previous section, is the coefficient of s\r{xi, . . . , Xk\zi, . . . , z n ) in the expansion 
of s\(xi, . . . , Xk\zi, . . . , z n ) x 1 (note the equality of second alphabets!). So it is equal to 5^ 
(and on the picture we have shown the unique A = A' configuration). 

The other regions are treated the same way as before. Regions D and E correspond 
to s\(xi, . . . , Xk\zx, . . . , z n ) and s M (xi, . . . , Xk\yi, • • • , y n )- If v encodes as before the edges 
between F and J, then region F corresponds to s v {x\, . . . ,Xk\yi, ■ ■ ■ ,y n ), on condition that 
Ai + fa < n — k. Finally region J is an upside-down MS-puzzle with the y^ in the reverse 
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Figure 14. Alternate graphical solution of the MS problem (n = 5, k — 3). 



order, so that we find the equality 

s A (xi, . . . ,x k \zx, . . . ,z n )s IJr (xx, . . . ...,y n ) 

= ^2 e A,^"' • • • > *i> • • • ' ^^(xi, . . . , Xfclyx, . . . , y n ) 

V 

As a corollary, we find the curious identity: 

z v \pky\i ■■■,y n ;z 1 ,...,z n )= e^(y n , . . . , yr, z u . . . , z n ) 

6.5. The Knutson— Tao problem. There is a more basic problem which consists in ex- 
panding the product of two factorial Schur functions with the same two sets of variables 
as a sum of factorial Schur functions, and more specifically providing a manifestly positive 
formulas for the structure coefficients in the sense of Graham. It was solved in various ways 
[13, 15, 10], but here we are particularly interested in the solution in terms of puzzles, as in 
[10]; this we call the Knutson-Tao (KT) problem. 

Clearly, a solution of the MS problem provides a solution of the KT problem by setting 
yi = Zi, i = l,...,n. Thus, the previous two sections provide two solutions of the KT 
problem, which turn out to be different. The most interesting one is the alternate one: 
when yi = Zi, the upper half of the rhombus J of the right hand side of Fig. 15 becomes an 
equivariant puzzle (as in the example shown), with the spectral parameters labelled in the 
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Figure 15. An example of tilings contributing to each side of the equality of Fig. 14. 
proper order, and the equality of section 6.4 becomes 

S\(xi , . . . , Xk | Z\ , . . . , Znjs^ {x\ , . . . , Xk | Z\ , . . . , Z n ^ ^ ^ ^ (zi , . . . , 2 n )S;y {x\ , . . . , Xfc | Zi , . . . , Znj 

v 

by definition of the c^\(zi, . . . , z n ), cf section 6.2. This is precisely the formula found in [10]. 

Appendix A. Square-triangle-rhombus tilings 

As mentioned in the introduction, it is known that the tiling model introduced in section 
2 is related to the square-triangle tiling model. More precisely, the latter is equivalent to 
the tiling model without the shaded tiles 7. We describe here the slightly more general 
(square-triangle- rhombus) tiling model which includes the shaded tiles (in the spirit of [2]). 
Though it is not needed anywhere in this paper, the connection is worth mentioning. 

Consider tiles of three types: equilateral triangles, squares and "thin rhombi" with angles 
30 and 150 degrees. All of them have sides of unit length. The square-triangle-rhombus tiling 
model consists in filling a region of the plane with these tiles, with an addition restriction on 
the allowed rotations of the thin rhombi which is the following. It is easy to see that all edges 
can occur in exactly six directions which differ from each other by multiples of 30 degrees. 
We may for example assume that they are of the form 15 + 30k degrees, k integer. Thus, 
there are four possible rotations of triangles, three rotations of square, and six rotations of 
thin rhombi among which we select only three as shown on Fig. 16. 
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a- a+ f3- (3+ f3 7_ 7+ 70 

Figure 16. Tiles of the square-triangle-rhombus tiling model. The colors are 
purely decorative. 




Figure 17. (a) The 6 directions of the generalized square-triangle model and 
how to deform them into the 3 directions (in blue) of the regular triangular 
lattice, (b) Example of square-triangle-rhombus tiling corresponding to an 
equivariant puzzle. 

Now consider the following operation on a given tiling configuration: deform it in such a 
way that the direction of every edge can only take three values, which is the closest slope 
of the form 60k, k integer. One can show that this can be done consistently (as always 
in the theory of random and quasi-periodic tilings, these tilings can be seen as projections 
of a two-dimensional surface embedded in a higher dimensional space - here, four - and 
the transformation corresponds to changing the direction of projection); furthermore, mark 
with a ± the edge which have been rotated ±15 degrees counterclockwise. In this way the 
triangles become one of the four tiles a of Fig. 1; similarly, squares become pairs of tiles of 
type /3; and thin rhombi become pairs of tiles 7. The operation is invertible and is thus the 
desired equivalence. As before, if we want to preserve integrability, we must allow only one 
of the three types of thin rhombi in any given region. 

For example, starting from Fig. 17 one obtains the equivariant puzzle of the right of 
Fig. 11. Note that one can read directly the three Young diagrams encoded in the three 
sides of the triangle (as was observed in [18] in the non-equivariant case), by completing the 
deformed puzzle into a hexagon with sides of lengths k,n — k, k,n — k, k,n — k and drawing 
the complement of the puzzle as sets of the three "forbidden" rotations of the thin rhombi. 
One must be careful that one recovers this way the three Young diagrams corresponding to 
the binary strings read clockwise. 

Another interesting observation is that the use of the Yang-Baxter equation (proposition 
1) corresponds in this new language to two elementary operations on tilings (cases (i,ii) and 
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(iii) in the proof), up to reflection and rotation: 




Note the strong similarity with the migration of [18]; the latter is roughly the special case 
of the Yang-Baxter equation where one spectral parameter is set to zero. 
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